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1. Introduction 

The Jordan canonical form (JCF) is a familiar canonical form under similarity of square matrices. 
It consists of a direct sum of Jordan blocks associated with eigenvalues, and it is unique up to per¬ 
mutation of these blocks [8, §3.1]. We assume throughout the paper that, for a given eigenvalue X, 
the Jordan blocks at X in the JCF are given in nonincreasing order of their sizes. In 1951 Flanders 
published the following result [4, Th. 2]: 

Theorem 1.1. If A e C mxn and B e C nxm , then the JCFs of AB and BA may differ only in the sizes of the Jordan 
blocks at 0. Moreover, the difference between two corresponding sizes is at most one. Conversely, if the JCFs of 
M e C mxm and N e C nxn satisfy these properties, then M — AB and N — BA, for some A, B. 

Theorem 1.1 has been revisited several times and re-proved using different techniques [1,10-12, 
14,16], In this paper, we investigate what happens if, instead of two matrices, we have products 

of k matrices, Ai, A 2 ,_ Ak e C nxn . We refer to products of A -[...., A/ ( , in any order and with no 

repetitions of the factors, as permuted products. 

We assume A\, _ Ak are all n xn to ensure all permuted products are well defined. An important 

difference between k = 2 and k > 2 is that, without any assumption on A]. Ak, the products of 

A!. Ak have, in general, completely different eigenvalues for different permutations. One exception 

is the eigenvalue 0: if 0 is an eigenvalue of some product of Ai,_A/ ( , then it must be an eigenvalue 

of every other product of A\,..., A/ { . Indeed, in Theorem 1.1 the eigenvalue 0 is treated exceptionally, 
with nontrivial results on the sizes of the Jordan blocks at 0. Flowever, the following simple example 
with k = 3 shows that the difference between the sizes of Jordan blocks at X = 0 can be arbitrarily 
large. 

Example 1.2. Let A = diag(l, 1/2,.... 1/n), B — — J n (—i) T , and C = (AB) -1 J n (0), where J n (X) is the 
nxn Jordan block at the eigenvalue X [8, Def. 3.1.1], Then ABC = J n ( 0) by construction, whereas CBA 
has a simple eigenvalue 0 and n — 1 nonzero eigenvalues. 

To verify the last statement, observe that B -1 is lower triangular with all elements on and below 
the main diagonal equal to 1. Therefore, the last two rows of B -1 A -1 are equal, up to the last-but-one 
entry. Flence, the last two rows of C = B -1 A -1 J n (0) are equal. Since CBA is a product of J n ( 0) 
with some nonsingular matrices, we have rankCBA = n — 1. Moreover, the vector Vo = [1,2,..., n] T 
belongs to its kernel because CBAvo = CB[1.1,..., 1] T = Cei = 0. Now suppose that there is a Jordan 
chain, so let V] be such that CBAv 1 = Vo. Then, since BA is invertible, there exists w such that 
Cw = Vo, but this is impossible as the elements of Vo are all different while the last two rows of C 
are identical. So by contradiction 0 must be a simple eigenvalue of CBA. 

Example 1.2 shows that it may be difficult to characterize the eigenvalues or Jordan block sizes for 
products of three or more matrices. 

In [3], Fiedler introduced a decomposition of an n x n companion matrix into a product of k = n 
matrices, C = n"=i M,-, and showed that the product of the matrices M, in any order is similar to C, 
hence all permuted products have the same JCF. For the nonzero eigenvalues, this is precisely what 
happens when k = 2, by Theorem 1.1. This motivates us to examine general conditions that allow an 
extension of Theorem 1.1 for nonzero eigenvalues to the case k > 2. The Fiedler factors M, have the 
following properties: 

(FI) Commutativity: MjMj = MjM if 1 1 — j| > 1. 

(F2) Mj are all nonsingular, possibly except for M„. 

Fiedler’s results suggest the possibility of extending Theorem 1.1 to products of three or more 
matrices under appropriate commutativity conditions. Indeed, we will show that if the graph of 
non-commutativity relations is a forest (see Section 4), then all permuted products have the same 
Jordan blocks for nonzero eigenvalues. This commutativity assumption generalizes condition (FI), and 
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imposes no requirement when k = 2, i.e., the two matrices can be arbitrary, thus recovering Theo¬ 
rem 1.1. We impose no nonsingularity condition such as (F2) because this imposes similarity, i.e., also 
the Jordan blocks at zero must be the same: an undesirable restriction given our goal of generalizing 
Flanders’ theorem. Indeed, Theorem 1.1 shows that the difference in the sizes of Jordan blocks at zero 
is at most 1 when k = 2. One key result of our paper is that, for general k, under our commutativity 
conditions this difference is bounded by k— 1, and the bound is attainable. 

For products of three matrices, our condition reduces to the requirement that one pair commutes, 
and we prove that the allowable sizes are exhaustive. More precisely, we prove that given two lists of 
these allowable sizes, there are matrices A, B, C such that the JCFs of ABC and CBA consist of Jordan 
blocks at X = 0 whose sizes match those in the respective lists. 

Several previous papers have addressed extensions of Flanders’ result to many matrices. For exam¬ 
ple, [7] examines cyclic permutations and [5] derives conditions for the products to have the same 
trace, the same characteristic polynomial, or the same JCF, with focus on k = 3 or 2 x 2 matrices. 
Unlike in previous studies, we deal more thoroughly with any permutation and arbitrary n and k ^ 3, 
and work with commutativity conditions guaranteeing that the Jordan structures for nonzero eigen¬ 
values coincide for all permutations. 

The paper is organized as follows. Section 2 reviews basic notions and previous results. In Section 3 
we analyze permuted products of k — 3 matrices. Section 4 discusses the case k > 3, which requires 
the use of permutations and graphs. We conclude in Section 5 with a summary and some open 
problems related to this work. 

2. Notation, definitions and some consequences of Flanders’ theorem 

We follow the standard notation /„ and 0„ to denote, respectively, the n x n identity and null ma¬ 
trices. Given a square matrix M e C nxn , A(M) denotes the spectrum (set of eigenvalues counting mul¬ 
tiplicity) of M; diag(/li,_ A/ ( ) is the block-diagonal matrix whose diagonal blocks are A\ . Ap, in 

this order (that is, the direct sum of At...., A^). Two matrices M, N e C nxn are similar if there is an 
invertible matrix P such that PMP -1 = N. 

The Jordan block of size k e N at zero is the k x k matrix 

TO 1 “I 


0 _ 

and the Jordan block of size katXeC is the k x k matrix Jp(X)Jp(0) + iU;<. 

For a given X e C, the Segre characteristic ofM at X, denoted by Sx(M), is the list of the sizes of the 
Jordan blocks at X in the JCF of M. In this paper we regard it as an infinite nonincreasing sequence 
of nonnegative integers, by attaching an infinite sequence of zeros at the end. Note that the Segre 
characteristic at any X is uniquely determined, and that this definition includes also those complex 
numbers that are not eigenvalues of M, though in this case all entries in the Segre characteristic are 
zeros. 

We use boldface for lists of nonnegative integers. Given two (possibly infinite) sequences of inte¬ 
gers ji = (jit , fi 2 , ■ ■ •) and fi' = (pt\, fi ' 2 ,...), we will often refer to the standard l°° and f 1 norms, 
which we denote by || • ||oo and || • ||i. 

Given k matrices At ...., Ap e C nxn , by a permuted product of AtAp we mean any of the 

products of At . A/ ( in all possible orders, without repetitions. The set of permuted products of 

A\, _ Ax is denoted by V(At, _A/ ( ). For instance, for three matrices A, B, C, we have 

P(A, B, C) = {ABC, ACB, BAC, BCA, CAB, CBA}. 

We will generally use the 77 symbol to denote elements of V(At .A/<). 

The following definition relates matrices M, N in Theorem 1.1, and plays a central role in this 
paper. 
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Definition 2.1. A pair of matrices (M, N), with M e C mxm and N e C nxn , is a Flanders pair if there are 
two matrices A e C mxn and B e C nxm such that M = AB and N = BA. In this case, we say that there 
is a Flanders bridge between M and N. 

We have the following elementary result: 

Lemma 2.2. IfM.Ne C nxn are similar, then (M. N) is a Flanders pair. 

Proof. If PMP -1 = N, with P nonsingular, then we may take B — PM, A — P -1 , which satisfy 
AB = M and BA = JV. □ 

The converse of Lemma 2.2 is not true in general. This is an immediate consequence of Theo¬ 
rem 1.1, since two matrices in a Flanders pair may have different Segre characteristic at zero and, as 
a consequence, different JCF. However, if M, N are nonsingular, then (M,N) is a Flanders pair if and 
only if M and N are similar. This is also an immediate consequence of Theorem 1.1. 

The relation TZ on C nxn x C nxn defined by “M1ZN if (M, N) is a Flanders pair” is not an equiva¬ 
lence relation, since TZ is not transitive. Moreover, Flanders pairs connecting three matrices M, N, Q_ 
in the form (M,N), (N, Q) are closely related to our problem. The following direct consequence of 
Theorem 1.1 establishes some elementary features of these pairs. 

Corollary 2.3. //Me (C mxm , N e C nxn and Q e C qxq are such that (M, N) and (N, Q) are Flanders pairs, 
then 

(i) Sx(M) =Si(Q),for all 0, and 

(ii) ||6b(M)-6b(Q)|| 00 <2. 

In Corollary 3.6 we give a characterization of pairs of matrices M, Q as in the statement of Corol¬ 
lary 2.3 and with the same size. We will see, in particular, that, when M and Q_ have the same size, 
the converse of Corollary 2.3 also holds. Corollary 2.3 can be extended directly to more than three 
matrices. 

Another feature of Theorem 1.1 we are interested in is its exhaustivity. The meaning of exhaustivity 
is exhibited in the following result. 

Theorem 2.4. Let p = (pt\, ptj, ■■ ■), and pt’ = (pt\, pt' 2 ,...) be two lists of integers with /ij 
and pt\ > ii' 2 ^ ^ 0, such that 

(i) ||/*.-jH/||oo<l, and 

(ii) Hplli =m, ||/t'|h = n. 

Then, there exist two matrices A e C mxn and B e C nxm , such that So(AB) = ft and Sq(BA) = 

Theorem 2.4 follows immediately from Theorem 1.1 just by noticing that it is always possible 
to construct two matrices M e C mxm and N e C nxn such that So(M) = fi and So(N) — p!, with 
m, n, p, p' as in the statement of Theorem 2.4. It can be proved also in a direct way by explicitly 
constructing A and B. This is the approach followed in [12, Th. 3.3]. We present an extension of 
Theorem 2.4 to three matrices in Theorem 3.4. Our approach owes very much to the one in 12], 

3. The case of three matrices 

Unlike what happens for two matrices, given three matrices, A,B,C eC nxn , the spectra of two 
different permuted products of A, B, C may be completely different. To verify this, one may just take 
three random matrices A, B, C and compute the eigenvalues of ABC and ACB. This is related to the 
fact that two similar matrices, as BC and CB are if one of B, C is nonsingular, may give two matrices 
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with completely different spectra when multiplied on the left by a third matrix A. To what extent 
may the spectra of different permutation products of three given matrices differ? One restriction 
is that the determinants must all be the same, which implies that if 0 is an eigenvalue of some 
permuted product then it must be shared by all permuted products. However, as we have seen in 
Example 1.2, the Jordan structure of the eigenvalue 0 may differ from one product to another. Let us 
first consider the case of nonsingular matrices. The following result shows that, without any additional 
assumptions, the only restriction on the spectra of permuted products of three nonsingular matrices 
A, B, C is that they all have the same determinant. It is a restatement, with a more straightforward 
proof, of Theorem 4 in [6j. 

Theorem 3.1. Let A i = {An,..., A„i) and A 2 — {Ai 2 , • • • > A n2 } be two sets of n nonzero complex numbers, 
with possible repetitions. IfX n ■ ■ • A„i = A 12 • ■ • X n2 , then there are three matrices A, B.C e C nxn , such that 
A(ABC) = A\ and A(ACB) = A 2 . 

Proof. By Lemma 2.2, it suffices to find two similar matrices M.N e C nxn , and a third matrix 
A 6 C nxn , such that A(AM) = 2 I 1 and A(AN) — A 2 . This can be done using only diagonal matri¬ 
ces. More precisely, set r\ ^0 (arbitrary), a\ :—X\\/r\ and, recursively for i = 2,..,,n, r,- \—Xi 2 la[-\, 
ai := An/r,-. Note that, with these definitions, we have 

(a 1 r 1 )(a 2 r 2 )---(a n r n ) X n X 2 i 

CZ n r 1 — — — A n 2 • 

(air 2 )(a 2 r 3 ) • • • (a n _ir n ) M 2 A 22 • • • ^n- 1,2 

Hence, if we set M = diagOr, r 2 , _ r n ), N = diag(r 2 , r 2 ,..., r n , ri), and A — diagfai. a n ), then M 

is similar to N, and AM — diag(An,..., A„i), AN = diag(Ai 2 . X n2 ), as required. □ 

Under the conditions of the statement of Theorem 3.1, by Theorem 1.1 we have A(ABC) = 
A(CAB) = A(BCA) = yli, and A(ACB) — A(BAC) — A(CBA) — A 2 . Moreover, as a consequence 
of Theorem 1.1, the set of permuted products is partitioned into two classes, namely: C\ — 
{ABC, BCA, CAB], and C 2 = { ACB. BAC, CBA). Any two products in each class are related by a “cyclic 
permutation”, so they form a Flanders pair. Hence, we can relate the JCFs of these permuted prod¬ 
ucts. The remaining question is to relate the JCFs of permuted products in Ci with the ones in C 2 . 
Theorem 3.1 shows that, if A, B, C are nonsingular, there may be no relationship at all between the 
spectra of products in different classes. 

Motivated by the work of Fiedler, here we require that at least two of A, B, C commute. As we 
see in Section 4, if we consider formal products of an arbitrary number of matrices, commutativity 
conditions allow us to characterize those cases where any two arbitrary permutations are linked by a 
sequence of Flanders bridges. In this case, all permuted products have the same Segre characteristic 
at each nonzero complex number. 

Proposition 3.2. Let A, B.C e C nxn be such that at least two of A, B, C commute. Let I7i, JJ 2 e V(A, B, C). 
Then 

(i) SxUL^^SiinriJorallX^O.and 

(ii) ||5 o(/7 1 ),5o(/7 2 )||o 0 <2. 

Proof. By Corollary 2.3, it suffices to show that, in the conditions of the statement, one of the follow¬ 
ing situations occurs: 

1. (/7i,/ 7 2 ) is a Flanders pair. 

2. There exists 77 e V{A, B. C) such that (IT \, IT) and (77, 77 2 ) are Flanders pairs. 

In the conditions of the statement there are at most 4 distinct elements in V(A , B, C), which give 
at most 6 distinct (non-ordered) pairs of permuted products. Let us assume, without loss of generality, 
that AC = CA. In this case, the elements in V(A, B. C) (including 77i and 77 2 ) are ABC, ACB, BAC, 
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CBA, and (ABC, ACB), (ABC, BAC), (ACB. BAC) and ( BAC,CBA ) are Flanders pairs. Hence, one of 
the situations described above holds for /7i and 772. □ 

The following technical Lemma 3.3 is used to prove Theorem 3.4: 

Lemma 3.3. Let fi — (/xi , /X 2 _), /x' = (pt\ , /x' 2 ,...) ef 1 be two sequences of nonnegative integers. Sup¬ 

pose that 

(i) ||/i — /t'lloo = 2, and 

(ii) ||/7||i = ll/7'Hr =n. 

Then we may rearrange /x and fi' in such a way that 


/7 (fli t , /7i 2 , /Xj 3 , /Xj 4 , /7i 5 , /Xj 6 , . ..), /7 (/Xq , /X( 2 , /Xj 3 , /X ( -^, /X^ , /X^, . . , 

with 


/7ij + /7i j+1 + /x ij+2 = n’i. + /x;. +i + n' jj+2 , for all j = 1 (mod 3). (3.1) 

Proof. Let m = max(||/7||o, ||/7'||o) be the maximum of the number of nonzero elements in /7 and /7 
We may assume that they both have the same length, by adding zeros to one of them if necessary. 
The proof is carried out by induction on m. For m ^ 3 the result is trivial. Suppose the result holds for 
lengths up to m — 1, and let us prove it for length equal to m. By condition (i) in the statement, there 
is some i > 0 such that |/x,- — /x'| = 2. Without loss of generality we may assume that /X; = /xi + 2. 
Now, condition (ii) in the statement implies that at least one of the following situations must occur: 

(Al) There is some j ^ 0 such that /x' = /Xj + 2; or 

(A2) There are some k, l ^ 0, with k ^ t, such that /xj ( = /x^ + 1 and \i! i = /x^ + 1. 

In case (Al), we may rearrange /7 and /x', by adding one extra zero in each list, in the form: 

/7 = (/Xj, /Xj, 0; /7), 

/7' = (/7<,/7'.,0; /x'), 

where /7 and jEt' are obtained from /7 and /x', respectively, by removing the ith and jth elements. 
Now, the result follows by the induction hypothesis on )it and /x'. 

In case (A2) we may rearrange: 


/t = (/xi,/x fc ,/xx;/t), 

/x' = (/x;,/Xfc,/X£;/t'), 

where /x and /x' are obtained from /x and /x', respectively, by removing the ith, kth and Fth elements. 
Again, the result follows by induction on /x and /x'. □ 

The main result of this section is an extension of [12, Th. 3.3] to three matrices A, B, C under the 
commutativity condition AC — CA. 

Theorem 3.4. Let fi — (/xi, /X 2 ,..., 0,...), /t' = (fi\ ,p! 2 , _0,...) ef 1 be two nonincreasing sequences 

of nonnegative integers such that 

(i) ||/t-/x'||oo ^2, and 

(ii) ||/x||i = ll/x'Hr = n. 
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Then, there exist three matrices A, B,C e C nxn , such that AC = CA and 
So{ABC) = fi, and So(CBA) = fi f . 

Proof. First, notice that if ||/t — /t/||oo ^ 1. then by Theorem 1.1 there exist A, B e C nxn such that 
Sq(AB) = ft and Sq(BA) = fi'. In this case, we may take C = I„ and we are done. Hence it remains 
to consider the case \\/i — fi'\\ 00 = 2. The proof reduces to showing that the statement is true in the 
following two cases: 

(Al) fi = (m, n, 0_), n' — (m — 2, n + 2,0,...), 

(A2) /l = (m,n,q, 0_), fl' — (m — 2,n + 1, q + 1,0_), 

with m, n, q ^ 0 and m ^ 2. Indeed, let us assume that the result is true for both (Al) and (A2), and 
let fi and \i! be as in the statement. By Lemma 3.3, we can rearrange fi and fi’ in such a way that 
they are partitioned as 

H = 0,...), and /i' = (n\,..., fi' a , 0,...), 

where the pairs (/£,,//,!) for i — are such that ||/t;||i = =: n, and they either sat¬ 

isfy ||fij — /ijlloo < 1 or are of one of the forms (Al), (A2). Now, since the result is true for 
both (Al) and (A2), and also for tuples of distance at most 1, there are matrices A\,B\,C\ e 
C" 1 xrli ,..., A a , B a , C a e C n “ xn “, such that A,C,- = C(A it and So(AiBjCi) = (fij, 0,...), S 0 (CiBiAi) = 
(/t', 0,...), for i = 1,_ a. Then the matrices 

A = diag(Ai,..., A a ), B = diag(£>i,.... B a ), C = diag(Ci,..., C a ) 
satisfy AC = CA and Sq(ABC) — ft, Sq(CBA) = 

It remains to prove that the result is true in cases (Al) and (A2). Consider (Al) first. Denote by 
Eq the matrix, of the appropriate size, whose (i, j) entry is equal to 1 and the remaining entries are 
zero. Set 

A = diag(/ m _i, 0, J n ), B = J m +n(0) + £ m +n,i > C = diag(0, 

Clearly we have AC = CA. Direct computation gives ABC = diag(/ m (0), J n ( 0)), and CBA — diag(0, 
Jm-2(0>, Jn+i(0)) + £ m +n,i. Now, diag(0, J n+ 1 (0)) + £„+ 2 ,i is similar to _/„+ 2 (0), because its only 
eigenvalue is 0 and its rank deficiency is one. Consequently, the JCF of CBA is diag(J m _2(0), J n+ 2(0)), 

so So(ABC) — (m, n. 0,...) and So(CBA) = (m — 2, n + 2. 0_), as required. 

Next consider (A2). Set 

A = diag(0, C = diag(/ n ^_q^_i, 0, / m _ 2 ), 

for which AC = CA, and set also 


B — diag(/q + i (0), Jm+n— l (0)) + £m+n+q, 1 • 

Direct computation gives 

ABC = diag(0, J q ( 0), J„( 0), J m -i(0)) + Em+n+q, 1, 


CBA = diag(jq +1 (0), / n+ i(0), J m _ 2 (0)). 

Note that diag(0, J q (0), J n (0), J m+ i(0)) + E m+n+qA is permutation similar to diag(_/ q (0), J n (0), 
diag(0, J m - 1 (0)) + £ m ,i). Since, as before, diag(0. J m _i(0)) + E m j is similar to J m ( 0), we conclude 
that So(ABC) = (m. n, q. 0,...) and So(CBA) = (m — 2, n + 1, q + 1,0,...), as required. □ 
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Remark 3.5. If ||/t — /t'||oo = 2, then the matrices A, B, C constructed in the proof of Theorem 3.4 
have the property that neither of the pairs (A, B) and (B, C) commutes, so there is exactly one 
commutativity relation in this case. In graph theoretical terminology (see Section 4), the graph of 
non-commutativity relations is a tree. 

Our last result in this section concerns the “non-transitivity” of Flanders pairs. 

Corollary 3.6. Let M, Q e C nxn . Then, the following conditions are equivalent: 

(a) There exists N e C nxn such that (M, N) and (N. Q) are Flanders pairs. 

(b) Sx(M) ~ S\(Q),for all XjLO, and ||So(M) - S 0 (Q)||oo «= 2. 

(c) There are three matrices A. B, C e C nxn such that AC = CA, M is similar to ABC, and Q is similar 
to CBA. 

Proof. The implication (a) =>■ (b) is Corollary 2.3. Suppose that (b) holds. Without loss of generality, 
we may assume that M and Q are given in JCF, so that M = diag(M r , M s ), and Q_ = diag(Q r , Q s ) 
where M r , Q r contain Jordan blocks associated with nonzero eigenvalues, and M s , Q s are Jordan 
blocks for X = 0. By hypothesis, we have M r = Q r and ||So(M s ) — So(Qs)lloo ^ 2. Using Theorem 3.4 
with fi = So(M s ) and /a' = So(Qs). we see that there exist A s , B s , C s such that A S C S = C S A S , 
A S B S C S = M s , and C S B S A S — Q s . The block-diagonal matrices A = diag(7 m ,A s ), B = diag(M r , B s ), 
C = diag(7 m , C s ), where m is the size of both M r and Q r , fulfill the conditions in (c). 

Finally, suppose that (c) holds. Let N = BCA. Then (M, N) is clearly a Flanders pair and, by the 
condition AC = CA, so is the pair (N, Q). □ 

We want to emphasize the difference between Corollary 3.6 and Theorem 1.1. The natural exten¬ 
sion of Theorem 1.1 to three matrices would be that (M, N) and (N, Q) are Flanders pairs if and only 
if there are three matrices A, B, C such that AC = CA and M = ABC, Q — CBA. However, we have 
not been able to prove that this is true and we have not found a counterexample. This issue remains 
an open problem (see Open Problem 3 in Section 5). 

4. More than three matrices 

For permutations in Ei ( , the symmetric group of (1.k}, we use the cyclic notation a — 

(iit 2 • •. i s ) to mean that er(ij) = i J+ i, for j = 1.s — 1, cr(i s ) — i i, and er(i) = i, for i ^ ij ,..., i s . 

An element in V(A\ . Ap) is related to a permutation a e Dp in the form A a - i ( i ) A 0 .-i( 2 ) • • • 

A„ t (k ), that is, er(i) is the position of the factor A, in the permuted product. In this case, we write 

n a '■= Ao-HvAa-iQ) • ’ ’ A a -i(ky 

Definition 4.1. Given a permutation a e Zk, a cyclic permutation of cr is a permutation of the form 
(1 2 ... k) e a, for some t ^ 0. We say that a, z are cyclically related if r is a cyclic permutation 
of a. 

Accordingly, given a permuted product TT a eV(A\ .A/ ( ), a cyclic permutation of IT a is a 

permuted product of the form 17 z e V(A\ . A\f), with r = (1 2 ... k) l o, for some l e N. 

If n a is a cyclic permutation of 77 T , then Tl a and 77 r are cyclically equivalent, and we write 

n 0 n z . 


We note that ~c is, indeed, an equivalence relation. Moreover, if /T cri ~c n a2 , then (77^,. n a2 ) is 

a Flanders pair. Conversely, if (/7 CTl , n a2 ) is a Flanders pair for all A-[ .A/ ( (that is, as a “formal 

product”), then 77^, is a cyclic permutation of 77^. 

Definition 4.2. Given two permutations oq, cr 2 £ 27 p, we say that i],... ,i g , with 1 < ij, ..., i g < k, 

have the same order in oq and 02 up to cyclic permutations if i 1 . i g appear in the same order in 

ff] := (1 2 ... k)“oq and 02 for some ct > 0. 
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Accordingly, given /7 CTl , n a2 e V{A\, _ Ap), we say that Aq,..., A; have the same cyclic order in 

both n„ x and n a2 if ii. i g have the same order in cri and (T 2 up to cyclic permutations. 

4.1. Inverse eigenvalue problem 

We start with an observation that characterizes Z^ up to cyclic permutations. 

Lemma 4.3. Let a, z e Z^ be two permutations. Then a and z are cyclically related if and only if all triples 
ii, 1 * 2 , i 3 , with 1 < ii, 12 , h < k have the same order in a and z up to cyclic permutations. 

Proof. If a — (1 2 ... k) e z, for some l ^ 0, then it is clear that each triple ij, i 2 , is has the same 
order up to cyclic permutations in both a and z. 

Conversely, assume that every triple ii, 12,13 has the same order up to cyclic permutations in a 
and z. Let a, f) > 0 be such that a := (1 2 ... k) a o and z := (1 2 ... k)^z satisfy <r(l) = 1 = r(l). 
Suppose a ^ r and let v = min{i: cr(i) ^ r(i)}. Then 1, a(v), z(v) do not have the same order up to 
cyclic permutations in a and r, a contradiction. Hence, a and r are cyclically related. □ 

We next show that it is possible that any two permuted products 77i, /T 2 have different spectra 
unless /7i ~c ^ 2 - 


Proposition 4.4. For each k ^ 3, there exist matrices, A 1 .A/ ( e C nxn such that for any two permuted 

products 77i and II 2 belonging to different equivalence classes of V(A\ . A\f) undergo, A(IT\) and A( 772 ) 

are different. 


Proof. First, let us order all the ( 3 ) triples ( 11 , 12 , 13 ), with 1 < ij < 12 < 13 ^ k using, for in¬ 
stance, the lexicographic order. This order induces an ordered list of length 3 ■ ( 3 ) = ^ ~ 2 ) , 

denoted by C, after adjoining all triples in the given order. For instance, for k — 4 we get the list: 
£ = (1,2.3; 1,2,4: 1,3,4: 2,3,4). Now, let y : {1, 2,..., W-iKfc~ 2) } (12 . fc} be the map de¬ 

fined by y (i) = C[ (the ith number in C). For each j = 1,.,., ( 3 ), by Theorem 3.1, there are three ma¬ 
trices Asj —2 , A3j_i, Asj e C 2x2 , such that A(Asj- 2 Asj-^Asj) i=- AiAs^As^-tAs^-s')- For i = 1. k, 

define 


At = diag(An, A i2 ,..., A ; q) e c 2 ^ 2 ©, 

where 


Aij = 


^3(j-l)+r, 

h. 


if there is some 1 ^ r ^ 3 such that y(3(j 
otherwise. 


l)+r) = i, 


For instance, for k = 4 we have Ai = diag(Ai, A 4 . A7, 1 2), A2 = diag(A2, A 5 , 12 , A10), A3 = diag(A3, /2, 
Ag, An), A4 = diag(/2. Ag, Ag, A12). 

Let n c r, and n a2 be two permuted products in V(A\,..., Aif) that are not cyclically equiva¬ 
lent. By Lemma 4.3, there is a triple ( 11 , 12 , 13 ), with 1 < iuh,h < k, such that i 1 , 12,13 appear 
in this order in a\, and they appear in the order 13 , 12,11 in as, up to cyclic permutations. The 

triple (ii, 12 , 13 ) corresponds to a triple (3j — 2,3j — 1.3j) in C for some j = l_ ,(3), such that 

A(Asj-2As]-\Asj) 7^ A(AsjAsj-tAsj-2). The result follows from the inclusions A(Asj-2Asj-\Asj) c 
A(n a ,) and A(AsjAsj-\Asj-2) ^ A(n a2 ). □ 


It is worth noting that, in the construction of the proof of Proposition 4.4, the spectra of /7 a , and 
n a2 are not necessarily disjoint. Note also that the size of the matrices, namely n = k(k — 1)(/c — 2)/2, 
depends on k. 

All permuted products in V(A\, ..., A/ ( ) have the same determinant. Equivalently, the product 
of their eigenvalues is the same for all permuted products. One may wonder whether this is the 





F. De Term et al. / Linear Algebra and its Applications 443 (2014) 120-138 


129 


only restriction on the eigenvalues of permuted products belonging to different classes under cyclic 
permutations, as it is for three matrices. More generally, we may pose the following problem. Here 
and hereafter, for a given set A of complex numbers, the notation X denotes the product of all 
numbers in A. 

Inverse eigenvalue problem for permuted products of k matrices: Given (k — 1)1 sets of n nonzero 
complex numbers, Ai,..., Ajk-iy, such that n^eyij ^ = 17 x 64 ^ A ' / or 1 ^ i> j ^ (k — 1)1, find matri¬ 
ces A], Ak, with Aj e C nxn , for i = 1 k, such that A(I7j) = Aj, for j = 1,.... (k — 1)!, where 

77j e V(A\,..., Ai<) belongs to the jth equivalence class under ~c- 

In Section 3 we saw that the “Inverse eigenvalue problem for permuted products of k = 3 matrices” 
is always solvable. The following result states that this is not true for k large enough. 

Theorem 4.5. Let n, k be two integers such that (k — l)!(n — 1) + 1 > kn 2 . Then, there exist (k — 1)1 sets of 
nonzero complex numbers A \...., A(k-iy, with | A,j = n and ^ — FLeyi K for all 1 0',j sj(k-l)!, 

such that there are no matrices A \,..., Ak e C nxn satisfying A(llj) = Aj, for j = 1.(k — 1)1, where 

77j e V(A \,..., A/ ( ) belongs to the jth equivalence class under Re¬ 
proof. We first note that prescribing the eigenvalues of a matrix A is equivalent to prescribing the 
coefficients of the characteristic polynomial paM := det(AJ — A). Let A\,..., Ak e C nx " be arbitrary 

matrices and let X = vec([Ai.A/<]) 6 C kn xl be the vectorization of the block matrix [Aj ... Ap] 

[9, Def. 4.2.9], Let us denote by TI\ . Ll(k-\)\ the representatives of each of the equivalence classes 

of V(A\,..., Ak) under ~c- Define the map 

P : C M -*■ C N , 

X (->■ P(X) = (Pj(X ),..., Pjv(X)), 

where P(X) is the vector containing the coefficients of the characteristic polynomials of Tl\ . 

n (k -\y, in a certain pre-fixed order. P is a polynomial map, since the coefficients of the characteristic 
polynomial of a matrix are polynomial functions of the entries of the matrix. Moreover, we have 
M = kn 2 and JV = (k — 1)1 (n — 1) + 1. Indeed, the necessary condition ri,eA, ^ ~ rixeyij ' A ' ^ or ^ ^ 
i, j ^ (fe — 1)1, is equivalent to the fact that the zero-degree coefficient of all characteristic polynomials 
of /7j, for j — 1,..., (k — 1)1, coincide. We may just slightly modify the definition of P, in such a way 
that, instead of n coefficients for each characteristic polynomial, we just have (n — 1) coefficients. 
Together with the choice of the determinant, this gives the (k — 1)1 (n — 1) + 1 coordinates in P(X). 

Now, the “Inverse eigenvalue problem for permuted products of k matrices with size n x n" is 
solvable, for k and n, if and only if P is surjective for these k and n. It is known that a polynomial 
map from C M to C N is not surjective when N > M [15, Th. 7, Ch. I, §6], so the result follows. □ 

4.2. Graph theoretical description ofV(A\ . Ak) 

We will see in Section 4.3 (see Theorem 4.16) how to characterize the maximum distance between 

the Segre characteristics of the zero eigenvalue in any two given products in V(A\ . A n ) using 

the graph of non-commutativity relations of pairs of matrices in {Ai. Ak}. From this graph, there 

arises an interesting combinatorial theory connected to this problem. The main result in this section 
is Theorem 4.13, which allows us to derive the main part of Theorem 4.16 as a direct consequence. 

For the basic notions in graph theory we follow [2], A graph is a pair of sets Q = (V, £), where 
V — (1,..., k} is the set of vertices, and £ is the set of edges', an edge is a two element subset of V. 
Here {i, j) e £ means that there is an edge joining i with j. By this definition, multiple edges between 
the same pair of vertices and edges joining a vertex to itself are disallowed. 

A sequence of edges {io, i'i}, {i'i, 12 },_ [i m - 1 . <m} is called a path of length m if all vertices ij are 

distinct. The sequence is called a cycle of length m if m ^ i m — i 0 , and all vertices ij, with 0 < j < m, 
are different from each other and i 0 . We say that a graph has a cycle if a subset of its edges is a cycle. 
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A graph Q = (V , £) is connected if, for any pair of vertices in V , there is at least one path containing 
them. A forest is a graph that has no cycles, and a tree is a connected forest. The degree of a vertex 
i e V in the graph Q — (V, £) is the number of vertices joined to i. A leaf is a vertex of degree one, 
and the parent of a leaf is the only vertex joined to it. A cut of G — (V, E) is a partition of V = Vi U V 2 
(Vi n V 2 — 0). Given a cut Vi, V 2 , we say that an edge [i. j} e E crosses the cut if i and j each lie in 
different Vp, b = 1,2. 

An oriented graph is a pair Q — (V, E) where V is again a set of vertices and £ is a set of ordered 
pairs of elements of V , that is, £ c V x V. Here (i, j) e E means that there is an edge joining i with 
j from i to j. A path (or cycle) of length m in an oriented graph is likewise a sequence of m elements 
of £ of the form (i 0 , i'i), (i 1 , 12 ), ■ ■ •, (i m - 1 , i m ). with distinct {ij} (except i m = i 0 for cycles). An acyclic 
oriented graph is an oriented graph with no cycles. 

An orientation of Q — (V , £) is a function co : £ -*■ V x V that assigns to each vertex [i, j } e E 
one of the ordered pairs (i, j) or ( j, i). Note that the set-wise image of £ under w, denoted a>(E) — 
(co(e): e e £}, associates a graph Q with an oriented graph (V,co(E)). The orientation to is said to be 
acyclic if (V,a>(E)) has no cycles. The set of acyclic orientations of Q is denoted by A(G). Any total 
order of V determines an acyclic orientation co e A(G) by &>({/, j}) = (i, j) if and only if i -< j and 
{/, j} e E. The converse is also true, as the following result shows. 

Proposition 4.6. Let Q — ( V, £) be a graph, and u> e A(G). Then, there is a total order =4ofV such that co is 
the acyclic orientation of Q determined by =<I (a topological sort of co). 

For a proof of Proposition 4.6 we refer the reader to [13, p. 137], 

To motivate the following definition, the graph that is our primary concern is Q = (V, £), where 
£ encodes the non-commutativity relations on k matrices (see Definition 4.14). An edge {/, j} e £ 
represents the fact that matrices A; and Aj do not commute. Meanwhile, we continue to associate 
elements of Xj. with elements of V(A\ . Aif) via n a . 

Definition 4.7. Given a graph Q — (V, £) with V = {1,2,_/<}, we say that to a is an allowed swap 

of a e X/ ( when r = (i i + 1) is a transposition with {cr —1 (t), cr -1 (i + 1)} £, for some 1 ^Uk-l. 

The proof of the following result is straightforward. 

Proposition 4.8. Let Q — (V, E) beagraphwith V = {1,2 ,,k}.Let~g be the relation on Dp defined by 

a\~q(J2 <£> 02 = T s o • • • o T2 o Tl o °T > 

where, for each i = 1,..., s, t,- is an allowed swap of T[_\ o ■ ■ ■ o q o o\. Then ~q is an equivalence relation. 

For a 6 Dp, we denote its equivalence class under ~g by [cr]g = {<? e Dp: a ~g a} and the 
quotient space (set of all equivalence classes) by Dp/ ~g. This set, as shown in Section 4.3, is closely 

related to the “generically” distinct elements in V(A\ __ Ap) required by the non-commutativity 

relations encoded in Q. From the combinatorial point of view, this set is interesting by itself because 
it is in one-to-one correspondence with the set of acyclic orientations of Q, as the following result 
shows. 

Theorem 4.9. Let Q = (\7, £) be a graph with V = {1,2,..., fc}. Let us define the map 
&G '■ Zk/ A(S), 

[er ]g T2g(cr), 

where, for each {i, jj e E, T2g(a) is the orientation given by 

&g(o 0({i, J}) := (i, j), ifcr(i) < cr(j). 

Then Qg is well defined (i.e., a ~g a implies £2g(cf) = Qg(cr)), and it is a bijection. 


(4.1) 
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Proof. Let us first show that fig is well defined. It suffices to show that it is well defined for a single 
allowed swap a — r o a, where r = (/ i + 1) and (cr _1 (i), er _1 (i + 1)} ^ E. Since r = (i i + 1), i] < h 
implies either r(ii) < r(i 2 ) or i\ — i 2 — 1 = i. Hence, cr(i\) < a( i 2 ) implies r o cr(ii) < r o cr(i 2 ) for 
all {ii, i 2 } e E since {o' -1 (i)» cr —1 (1 + 1)} ^ E. 

Now, let us show that fig is surjective. Let to e A(Q) be an orientation of Q. According to 
Proposition 4.6, there is a topological sort that produces a total order =<: on V, which we write as 

h -< h -< • • • -< ik- Then we take a e Zp defined by cr(j ) = ij, for j = 1 . k, so that to = fig(cr). 

To prove that fig is injective, let a and a be two permutations in X), such that fig(a) — fig (a). 
Without loss of generality, and relabeling the vertices of V if necessary, we assume that a = id is the 
identity permutation. We construct a sequence of permutations, a\ — a, er 2 .ff/<, by the recurrence 

of+i = {i i + 1 ... cr,-(i)) ocr;, (4.2) 

for 1 < i < k (note that this recurrence implicitly requires showing er,(i) ^ 0- The proof reduces to the 
following two claims: 

(i) err 1 (j) = j, for j < i (in particular up = id), and 

(ii) o - i+1 ~g CTi. 


Note that the first claim also implies both cr,-(j) ^ i and cr~^(j) ^ i when j ^ i, which justifies the 
requirement that a,(i) ^ i. 

We proceed by induction on i. For i = 1, both claims are trivial. We now assume both (i) and (ii) 
are true up to some i < k. If it happens that er, (D = i, then both claims are trivially satisfied at i + 1 
with er; + i = cr,-. Otherwise, for cr,_,_i as in (4.2), we have, for the first claim: 

for j <i: (r~ + \(j)=a~ 1 o (03 (i) ... i)(j) 

= cr-\j) = j, 

for j = i: a~\ (i) = a. 1 o (ct,( i) ... i ) (i) 

= txf 1 (or/CO) = i. 

For the second claim, let cr,-j = (j ... crKO) 0 cr/, for i j < cr,-(i). Note that crrA 1 (j) = cr~^(j) and 
ar'(j + 1) = cr-\orj(i)) = i. Then, ary = (j j + 1) o< 7 y +1 is an allowed swap unless (rrr 1 (j), i} e E. 
But we have that — j < er,(i) and, on the other hand, by claim (i), we have cr~^(j) > i. 

Hence {rr^ 1 (j), i} ^ £, since fig(<7t) = ^2c;(ct) = fiq(id), by the induction and the initial hypothe¬ 
ses. □ 

According to Definition 4.1, we introduce the following notion. 

Definition 4.10. Given two classes Ci and C 2 in Zp/ ~g, we say that Ci and C 2 are cyclically related if 
there are some cr, e C\ and 02 e C 2 such that o\ and 02 are cyclically related. 

We note that, unlike the relation for permutations in Zk introduced in Definition 4.1, the relation 
on equivalence classes in Definition 4.10 is not necessarily transitive. 

We show in Theorem 4.12 that the cyclic relation of a pair of elements of Z]J ~g corresponds to 
the following relation between the corresponding elements of A(Q). 

Definition 4.11. Let Q = (V, E) be a graph and A(G) be the set of acyclic orientations of Q. Given 
£ A(Q), we say that to 2 is a cut-flip of if there is a cut V = V] U V 2 such that, for each 
edge {i, j} e E, we have: 

(a) &>i({i. j}) = co 2 ({i, j}) if i, j e Vi or if i,jeV 2 . 

(b) cordi, j}) = (i, j ) and to 2 ({l j}) = U, 0 if i e Vi and j e V 2 . 
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We say that co,coe A(G) are connected by d cut-flips if there exists a sequence — a>, a >\,..., — 

w e A(G), such that &>,-_! is a cut-flip of &>,, for i — 1. d. 

We say that A(G) is connected by d cut-flips if any two orientations are connected by d cut-flips. 
For d = 1 we just say that A(G) is connected by cut-flips. 


We note that, by swapping Vi and V 2 , the relations in Definition 4.11 are symmetric. In plain 
language, a cut-flip is a cut where the edges of Q that cross the cut are oriented from V\ to V 2 in 
and from V 2 to V\ in a >2 , while the non-crossing edges of Q have the same orientation in both w\ 
and &> 2 - Theorem 4.12 shows that cut-flips graphically represent the cyclic relations of quotient space 

Zk/~G- 


Theorem 4.12. Let Q = (V, E) be a graph, with V = {1.2,...,/<}, and let Qg be the map defined in Theo¬ 
rem 4.9. Then [oq]g and [oqlq ore cyclically related if and only if T2q ( 02 ) is a cut-flip of Qq(t Ti). 


Proof. For brevity, throughout the proof we set r := (1 2 ... k) e . For a e we have 


rod (i) = 


cr(i) + l, if cr(i)<k — l, 

cr(i) + t — k , if cr(i)>k — £, 


and hence 


r o er(j) — r o cr(i) 


a(j) - a (i), 
< 0 , 

> 0 , 

cr(j) - cr(i), 


if cr(i), a(j) ^ k — i, 
if er(i) ^ k — £ < cr(j), 
if er(j) ^ k — i < cr(i ), 
if o'O'), a(j ) > k — l. 


(4.3) 


Now, let us prove the “only if’ part of the statement. Let a 6 [oq ]q and a e [a 2 \q be such that 
a — t o a, for some l. Recall that, by Theorem 4.9, Qq(cr) — Z2q(o\) and Qqfa) = 12q(02). Let us 
consider V\ = {er -1 (l),.... cr _1 (k — £)} and V 2 = {cr —1 (k — l + 1 ),..., er -1 (k)}. 

We now verify that Qq(a) is a cut-flip of Qq(cr). We analyze all possible situations for an edge 

0J}e£: 


• i, j e V]: Then, o(i), a(j) ^k — l, by the definition of V\. Hence, by (4.3), o(j) — o(i) = a(f ) — 
o'O’), so T 2 g(a)({i, j}) = T 2 g(a)({i, j}). 

• U j £ V 2 '. Then, a(i), a(j ) > k — i, by the definition of V\. Again, by (4.3), a(j) — a(i) = a(j) — 
o(i), so S2g(0)({i, j}) = S2g(cr)({i, j }) as well. 

• i e V\, j e V 21 Then cr(i) ^ k — l < o'O'), by the definition of VO and V 2 , so cr(i) < cr(j). Also, by 
(4.3), a(j) -<r(i) < 0, so S2g(a)({i, j}) = (j, i), whereas S2g(cr)({i, j}) = (i,j). 

• i e V^ 2 , j G V\'. In this case, er(j) ^ k — l < cr(i), so a{j) < cr(i). Again, by (4.3), <f(j) — cr(i) > 0, 
so %(ff)({i, j}) = (i, j), whereas Qg(cr)({i, j}) = ( j , i). 


Let us now prove the “if’ part. Suppose T2g(o\) is a cut-flip of £2g(o 2 ) via the cut Vt U V 2 . Set 
l k — |V0|. Let 0 be the total order of V defined by: i =4 j if and only if aqO') ^ oq(j). Define 
a permutation tt e Zk by having tt - 1 (1), ..., 7t _1 (k — T) be the elements of Vi sorted according 
to and jt _1 (k — £ + 1),..., 7r _1 (k) be the elements of V 2 also sorted according to =<;. Note that 
7 t(i) ^ k — l < 7 t(j) for all i e Vj, j e V 2 . 

From Theorem 4.9, the problem now reduces to showing 

T2g(n) = T2g(a-[) and £2g(x on) = &g( 02 ), 

or equivalently, for all [i, j} e E, n(i) < n(j) if and only if oq(i) < oq(j), and r o n(i) < r o n(j) if 
and only if oqO) < a 2 (j)- As before, we consider the separate cases: 


• i, j e : Then n(i), n(j ) ^ k — l, and, by the definition of n, n(i) < n(j) if and only if i -< j, 
which in turn holds if and only if oq(i) < oq(j), by the definition of =<:. Also, by (4.3), r o n(j) — 
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t o 71 (i) = 7 r(j) — 7r(i). Note that, by the definition of cut-flip, the sign of this difference is in turn 
equal to cr 2 (j) - cr 2 (0- 

• i, j eV 2 '. Similar arguments lead to 7 t{i) < 7 T(j) if and only if eri(i) < cri(j) and to jz(i) < to 7 t{j) 
if and only if 02 (0 < 02 (j) also in this case. 

• i e V\, j £ V 2 : In this case, we have 7 t(i) < 7 z(j) by construction and oq (i) < oq(j) by hypothesis. 
Also, r o 7 t(i) > z o 7 T(j) by (4.3) and, by the definition of cut-flip, cr 2 (i ) > cr 2 (j). 

• i e V 2 , j e Vi: With similar arguments, we have tt(i) > 7 t(j), 01 ( 1 ) > 0 i(j), and r o jr(i) < z o 
7t(i). 02 ( 1 ) < 02 0')- □ 

The main result in this section is Theorem 4.13, which gives a simple characterization of those 
graphs Q such that Ek/ ~g is connected by cut-flips, and, when this is the case, establishes the 
maximum number of cut-flips needed to connect any two classes. 

Theorem 4.13. LetQ = (V, E) be a graph with k vertices. Then A(G) is connected by cut-flips if and only if Q 
is a forest. Furthermore, ifQ is a forest, and d < k is the length of the longest path in Q then any two classes in 
A(Q) can be connected by no more than d cut-flips. 

Proof. Proceeding by contradiction, we assume that A(Q) is connected by cut-flips and that there 
is a cycle in Q, given by {io, ii), {ii, i 2 ],..., {i m _i, io) e E. Let a>i and co 2 be two acyclic orientations 
of Q related by a cut-flip. Every cut must be crossed by an even number of edges in the cycle, so 
the number of edges of the cycle on which a>\ and co 2 disagree (that is, the number of j for which 
ij}) ^tt> 2 ({ij_i. ij})) must be even. Hence, any two acyclic orientations with an odd number 
of disagreeing edges on the cycle cannot be connected by a sequence of cut-flips. Since there always 
exist two such acyclic orientations, we get a contradiction. 

To prove the converse implication, let us assume that Q is a forest. We will prove, by induction 
on d (the length of the longest path in Q), that Q is connected by, at most, d cut-flips. 

For d = 0, A(Q) has only one element, since no edges means that there is only one (vacuous) 
orientation. For d= 1, given two different acyclic orientations co,co of Q, we take Vi to be the set 
of all i where {i, jj e E for some j and co([i, j}) — (i, j) co([i, j}) — ( j, i) with V 2 = V — Vi. By this 
construction, all edges in E where co and w agree join vertices which are both in V 2 and each edge 
where they disagree is oriented from V 1 to V 2 by co and from V 2 to V\ by S. 

We now assume the result for d. Let L be the set of leaf nodes of Q. Consider the graph Q = 
(V — L, E) obtained from Q by removing L from V with a subset of the edges E = {{/, j} e E: i. j e 
V — Lj. Q is a forest with longest path length d — 2 (since any maximal path in Q must start and end 
on leaf nodes). Let co and S be two different acyclic orientations of Q. Then w and co, as functions 
restricted to E, are orientations of Q. By induction, co |-g and co are connected by at most d — 2 
cut-flips, So = co Ij, S i__ cbq = S|j for q ^ d — 2. For each S p in the sequence, we define an orien¬ 

tation co p of Q by extending S; to a function on E taking co p | £ _j = co | £ _j (in other words, co p agrees 
with co on the edges of the leaf nodes). As Sj_i,S, are related by a cut-flip with cut V\, V 2 (where 
Vi U V 2 = V — L), C 0 j-i, cot are related by a cut-flip with cut = V] U Li, V 2 — V 2 U L 2 where Li is 
the set of all leaf nodes whose parents are in V\ and L 2 = L — L\. 

Hence, co q is connected to co by q < d — 2 cut-flips, and co q \ p = co\ p , co q \ E _ P = co\ E _ p . What remains 
is to connect co q to w with 2 cut-flips. We identify two disjoint subsets of L corresponding to whether 
coq orients the leaf first or second in its edge: 

Mi = {ieL: {i,j}eE, co q ({i, jj) = (i, j) /S({i, j}) = (j,i)j, 

M 2 = {i e L: {i, j] e E, co q ({i, j j) = (j, i) / S({i, j }) = (i, j)}. 

The cut-flip given by Mi, V — Mi followed by the one given by M 2 , V — M 2 , connects co q to co. □ 

By Theorem 4.9, via the map Qg we may identify the quotient space EiJ ~g with the set of 
acyclic orientations of Q. Theorem 4.12 tells us that A(Q) is connected by cut-flips if and only if any 
two equivalence classes in Ek/ ~g are cyclically related. As a consequence, Theorem 4.13 says that 
any two classes in E]J ~g are cyclically related if and only if Q is a forest. 
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4.3. Commutativity conditions and distance of Segre characteristics 

As shown in Proposition 4.4, when there is more than one equivalence class in V(A\ .A/<) 

under ~c, it is pointless to ask about the change in the JCF of different permuted products, since the 
spectrum can be completely different. On the other hand, when there is only one equivalence class, 
all permuted products have the same nonzero eigenvalues with the same Segre characteristic. 

Motivated by Fiedler matrices, we impose certain commutativity conditions such that any two 

elements of V(A\ .A; ( ) are connected by a sequence of Flanders bridges, as we did in Section 3 for 

three matrices. Under these conditions, we also analyze the change in the Segre characteristic of the 

eigenvalue zero for different permuted products. We say that two products 77 and 77 in V(A\ .A;<) 

are related by a sequence of Flanders bridges if there are some ri\ ./7d+i e V(A\ . Ak) such that 

77i = 77, /7d + i = 77 and (77,, 77, + i) is a Flanders pair, for i — 1. d. 

Our commutativity conditions are encapsulated in the associated graph. More precisely, we are 
interested in the graph comprising the non-commutativity relations. 

Definition 4.14. Given k matrices A\,...,Ak e C nxn , the graph of non-commutativity relations of 

A\ __ Ak is the graph Q — (V, E) with V = {1.2such that {i, j } e E if and only if AjAj ^ 

AjAi, for all 1 < i, j < k with i ^ j. 

Given k matrices A\,...,Ak e C nxn , the set of permuted products V(A\ .A/U can be ana¬ 

lyzed in the light of the combinatorial approach of Section 4.2. In particular, if Q = (V , £) is the 
graph of non-commutativity relations of A\, _A/ { , and a e then an allowed swap of a ex¬ 
changes two consecutive factors in TT a , h> A a i(i + i ) A rJ i (j y This is allowed because 

{cr _1 (i), er _1 (i + 1)} ^ E, which means that the matrices A a - i (i ) and A a commute. Hence, the 

equivalence classes in the quotient space Epf ~g correspond to the generically distinct elements 
of V(A \,..., Ak) obtained by the commutativity relations of the complementary graph of Q. Here 
the word “generic" means that, for some particular A\,...,Ak, it may happen that some permuted 
products in V{A \,..., A/<) coincide even if they belong to different equivalence classes. For instance, 
ABC = ACB is possible even if BC ^ CB, though this is not generically the case. 

Hence, if we consider the elements of V(A\ .A/ ( ) as “formal products”, that is, as words of k 

letters, Aj,..., Ap, then [criJg and \G 2 \q are cyclically related if and only if there is a Flanders bridge 
between 77 CTl and n a2 . Then, using Theorem 4.12, the first part of Theorem 4.13 can be stated as 
follows. 


Theorem 4.15. Let A\, A/< e C nxn and Q be the graph of non-commutativity relations of A] . A/ { . 

Let permuted products in V(A\,..., A/U be considered as formal products (that is, as words of k letters, 

Ai,_ Ak). Then any two products in "P(Ai.A/ ( ) are related by a sequence of Flanders bridges if and 

only ifQ is a forest. 

Now, from Theorem 4.15 and part (i) of Corollary 2.3, we conclude that when Q is a forest, all 
permuted products in V(A\ .A^) have the same nonzero eigenvalues together with their cor¬ 

responding Segre characteristics. The remaining question is to analyze what happens to the zero 
eigenvalue in this case. Theorem 4.16, which is the main result in this section, establishes an up¬ 
per bound for the distance of the Segre characteristic at zero of two permuted products, and shows 
that the bound is attainable. This bound comes from the number of Flanders bridges in the sequence 
that relates two arbitrary permuted products 77^ and n a2 ; that is, the number of cut-flips connecting 
Qg(cri) and Qg(a 2 ). 


Theorem 4.16. Let A\, A/, e C nxn and Q be the graph of non-commutativity relations of A] .A/ { . 

Assume that Q is a forest and let d be the length of the longest path in Q. Then, given 77i, 772 e V{A\ ,..., A/ ( ), 
we have 


5o(/7 1 )-5 0 (77 2 )| oo <d. 


(4.4) 
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Moreover, this bound is attainable in the following sense: Let Q be any forest with k vertices, and let d < k 
be the length of the longest path in Q. Then there exist k matrices At,..., Ak such that Q is the graph of 
non-commutativity relations of At,..Ak, and there are Tit, Tlje V(A\ . Ak) with 

!<So(77 1 )-SoU7 2 )|| 00 = d. (4.5) 

Proof. The first part of the statement is an immediate consequence of Theorem 4.13. More precisely, 
let Q, TI\, and 77 2 be as in the statement. By Theorem 4.13, A(G) is connected by at most d cut- 

flips. This implies, using Theorem 4.12, that there are < 71,02 . a d+i e such that 77i = 77^,, 

TTa d+1 = 77 2 , and such that (TT aj , n aj+1 ) are Flanders pairs, for i = 1,_ d. Now, Theorem 1.1 gives 

|| S 0 (n ai ) - So(77 CTj+1 ) | ^ ^ 1. for i = 1. d, 

so we have 

d 

|| So (Hr) - So (77 2 ) || „, < £ || S 0 (/7 ffj ) - S 0 (/7 fff+1 ) || ^ < d. 
i =1 

For the second part of the statement (attainability of the bound (4.4)), we first consider the case 
where Q is a path of length d, and define the (d + 1) x (d + 1) matrices 

Aj = diag(/ d _i, J 2 ( 0), Ij-t), for i = 1,..., d, 

A d +t = diag(0, I d ). (4.6) 

The graph of non-commutativity relations of A\, _A d+ i is a single path of length d from A\ 

to Ad+t- Moreover, we have 77i = At A 2 ... A d+ i = y d+ i(0), and 77 2 = A d+ i... A 2 Ai = 0 d+ i, so 
l|So(7Ti) — S 0 (/7 2 )||oo = d. 

If Q — (V, E) is a tree with V = {l,...,k}, let us assume, without loss of generality, that 

{1,2), {2,3),..., [d, d +1} is a path of length d in Q. Now, let At, _ A d+ i be as in (4.6), and A d+2 = 

• • • = Ak = l d+ t- Let us number the edges in Q different from (1,2}. [2, 3}.(d, d-l- 1}, as et . e g . 

For each of these edges we build up the following k matrices: for the edge e s = {i, j), with 1 =% s < g, 
let D^\ ..., D[ s) be k nonsingular matrices of size 2x2 such that d[ s, D® ^ Dj S) Dj s \ and Dg S> = 1 2 
for l ^ i, j. Now, set At = diag(Aj, Dj 11 ,_ d| s> ), for i — 1. k. The graph of non-commutativity re¬ 
lations of Ai,_A/ ( is Q, by construction. Moreover, since D| s) is nonsingular, for all s = 1. g and 

i = 1 — ,k, we have 77! :=AtA 2 ...A k = diagf _/ d+1 (0). Mi) and 77 2 := A k ... A 2 Ai = diag(0 d+ i, M 2 ), 
with Mi, M 2 nonsingular, so So(TTt) — (d +1), and So(77 2 ) = (1,..., 1) (containing d +1 ones), hence 
||5o(77i)-5o(/7 2 )|| 00 = d. 

Finally, let G be a forest with t trees. Let l<t ,..., k t be the number of vertices in each tree, 

with kt + • • • + k t — k, and let dt, _ d t be the lengths of the longest path in each tree. By 

hypothesis, we have max{dj: j = 1 t} = d. For each tree, say the jth one, we define matri¬ 
ces Aj J) A^ ] 6 C"J xn ' as before, such that the graph of non-commutativity of A^ J) A^ 

is precisely this tree, and such that ||<So(A^ J> A^’... Aj/ 1 ) — SofA^ ... A^A^'lHoo = dj. Now, we 

set A, = diaglA^ 0 .A^), for i = 1. k, where A ( . l} = A^\ if i = kt + ■■■ + fcj-i + h, for 

some 1 ^ h ^ kj (where we set ko := 0), and A^’ = I nj otherwise. For these matrices, we have 

l|S 0 (A, A 2 ■ ■ ■ A k ) - So(A k ■ • ■ AzAijUoo = max , ,.. \\So(a\ ]) A ( 2 j) ... A^) -<S 0 (A^ ... A‘ J) a“ , )|| 00 = d, 

l|5o(Ai A 2 ... Ak) — So(A/ ; ... A 2 Ai)||oo = max J= i. t dj = d, and the graph of non-commutativity rela¬ 

tions of At,..., Ak is G, by construction. □ 

The construction in the proof of Theorem 4.16 does not necessarily give the minimum size of 

Ai. Ap that satisfy the second part of the statement. Also note that d^k— 1, and equality holds 

if and only if G is a path of length k — 1. 
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We mention that, in the case in the Fiedler matrices [3] the graph of non- 

commutativity relations is a forest. Moreover, it is just a path from Mi to M„. 


Example 4.17. Let Q — (V, E), with V = {1,2.9}, be the following graph: 



The length of the longest path in Q is d = 4, which corresponds, for instance, to the path 
{9,1}, {1,3}, {3, 8 }, { 8 , 7}. 

Now, let us follow the construction in the second part of the proof of Theorem 4.16 to 
get 9 particular matrices A\,...,Ag such that Q is the graph of non-commutativity relations of 
Ai,...,Aq (by identifying the jth vertex of Q with the matrix Aj) and so that the products Tl\ = 
(AqA\AiA^Ai)AqA 2 AsA^ and II 2 = (AtA^A^AiA^A^AjA^Ah satisfy ||<So(7Ti) — <So(^ 2 )lloo =4. 

Set A 9 = diag(/ 3 , ] 2 ( 0)), A\ = diag(/ 2 , J 2 ( 0 ), 1), A 3 = diag(l, J 2 ( 0), 1 2 ), A s = diag(J 2 (0), 1 3 ), A 7 = 
diag(0, 14 ), and A, = I 5 , for i ^ 1,3,7, 8 , 9. Now, let us number (and label) the edges that are not 
in the path {9,1}, {1,3}, {3, 8 ), { 8 , 7} as follows: ej ={3,6}, e 2 = {2.3}, e 3 = {3,5}, e 4 = {4, 5}, and 
set 


Ai = diag(Aj, / 8 ), A 2 = diag(/ 7 , of, / 4 ), A 3 = diag(A 3 , d' 1 ’, D< 2) , D* 3) , J 2 ), 

A 4 = diag(/ 11 ,D' 4) ), A 5 = diag(/ 9 , Df \ Dg 4) ), A 6 = diag(/ 5 , d' 1 ’,/ 6 ), 

A 7 =diag(A 7 ,D‘ 2) ,/ 4 ), A 8 = diag(A 8 ,J 8 ), A 9 = (A 9 ,/ 8 ), 

with D ( ‘ ] being nonsingular 2x2 matrices such that D^Dg 11 ^ Dg^Dg 11 , D^D^, 

D 3 3 ) D ( 5 3> ^ DfDf, and D l 4 ) D <4> ^ D ( 5 4 > D ( 4 4> . Under these conditions, the graph of non-commutativity 
relations of A\, _Ag is Q, and we have 

ri\ = {AgA\AiA%A 2 )AsA 2 A 5 A 4 = diag(js( 0 ), /), 

and 


n 2 = (A-jAp,A-iA\Ag)A^A 2 A 3 A 4 = diag(0s, J), 

with J = diaglD^’Dg 1 ', D^D®, D 3 3 > D ,3 \ D i 4 ) D >41 ). Since the matrices D 1 . 11 are nonsingular, we have 
that «S 0 (/7i) = (5) and «S 0 (77 2 ) = (1,1,1,1,1), so ||5 0 (/Ti) - 5 0 (/7 2 )||oo = 4. 

Theorem 4.13 implies that if the graph of non-commutativity relations of A\ . Ai ( has no 

cycles, then all permuted products in V(A\ . Ap) have the same eigenvalues with their cor¬ 

responding Segre characteristics. The reverse implication, however, is not true. For example, take 
A],..., Ape C nxn to be upper triangular and such that no pair commutes and the products of the 
( 1 , 1 ) diagonal entries of all matrices, 7Tj = Ai(i, i)A 2 (i, i)... Ap(i, i), satisfy jTi^TTj for i ^ j. Then, 
all permuted products have the same eigenvalues, with the same Segre characteristic (they are all 
simple eigenvalues). However, the graph of non-commutativity relations is the complete graph with k 
vertices, which is far from a forest. 

5. Conclusions and open problems 

In this paper we have analyzed the change in the JCF of products of k square matrices under per¬ 
mutations of the factors. As an immediate consequence of a classical result by Flanders, the products 
are classified into equivalence classes under cyclic permutations of the factors, in such a way that 
the structure in the JCF for nonzero eigenvalues coincides in any two products belonging to the same 
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class. We have first shown that, if no assumptions are imposed on the factors, then any two prod¬ 
ucts belonging to different classes under cyclic permutations may have different nonzero eigenvalues. 
Moreover, for three matrices, we have seen that it is always possible to prescribe the nonzero eigen¬ 
values of ABC and ACB, with the condition that the product of all eigenvalues coincides for both 
products. However, we have seen that this prescription is not always possible for more than three 
matrices. 

We have further shown that, by imposing certain commutativity conditions on the factors, the 
structure in the JCF corresponding to nonzero eigenvalues coincides for all products. In particular, 
we have seen that this is always the case if the graph of non-commutativity relations of the factors 
is a forest. We proved that, when considering formal products, there is only one equivalence class, 
up to cyclic permutations of the factors, if and only if the graph of non-commutativity relations of 
the factors is a forest. When this graph is a forest, we obtained an upper bound on the difference 
between the structure (sizes of Jordan blocks) of the JCF associated with the eigenvalue zero in any 
pair of products, and we saw that this bound is attainable. Moreover, in the case of three matrices, 
we proved that it is always possible to prescribe the sizes of the blocks associated with zero in the 
JCF of ABC and ACB as long as the difference between the corresponding sizes is at most 2. 

We conclude with some open problems that arise as a natural continuation of the problems ad¬ 
dressed in this paper. 

• Open problem 1: Is it always possible to prescribe the n eigenvalues of the (k — 1)! classes under 
cyclic permutations, provided that the product of all eigenvalues is the same, for k, n satisfying 
(k — l)!(n — 1) + 1 ^ fcn 2 and k > 4? 

• Open problem 2: Given d ^ 4 and two nondecreasing sequences //,, fi' of nonnegative inte¬ 
gers such that ||fi — /t'Hoo ^ d — 1, is it always possible to find d matrices, A\, _A^, such 

that their graph of non-commutativity relations is a path, and such that SofA!... Ad) = /*, and 
<So(Arf ... Ai) = /*,'? (The extension of Theorem 3.4 to d ^ 4.) 

• Open problem 3: If JW, Q. e C nxn are such that S\(M) = Sx(Q), for all X ^ 0, and ||5o(JVf) — 
S 0 (Q)||oo SJ 2, are there three matrices A, B,C e C nxn with AC = CA, such that M = ABC and 
Q = CBA? 

• Open problem 4: Obtain necessary and sufficient conditions for all products of a given set of k 

matrices to have the same nonzero eigenvalues and corresponding Segre characteristics (in the 
notation of the paper: Sx(Tl^) = Sx(f[ 2 ), for all X =A 0, and for all /7 i,/72 eV(A\ .A/ ( )). 
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